Enhanced Sensitivity to the Time Variation of the Fine-Structure Constant and nip/ me 
in Diatomic Molecules: A Closer Examination of Silicon Monobromide 
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Recently it was pointed out that transition frequencies in certain diatomic molecules have an 
enhanced sensitivity to variations in the fine-structure constant a and the proton-to-electron mass 
ratio mp/m^ due to a near cancellation between the fine-structure and vibrational interval in a 
ground electronic multiplet [V. V. Flambaum and M. G. Kozlov, Phys. Rev. Lett. 99, 150801 
(2007)]. One such molecule possessing this favorable quality is silicon monobromide. Here we 
take a closer examination of SiBr as a candidate for detecting variations in a and m^jm^. We 
analyze the rovibronic spectrum by employing the most accurate experimental data available in the 
literature and perform ah initio calculations to determine the precise dependence of the spectrum 
on variations in a. Furthermore, we calculate the natural linewidths of the rovibronic levels, which 
place a fundamental limit on the accuracy to which variations may be determined. 



PACS numbers: 06.20.Jr,33.20.Bx,33.15.Mt,33.15.Pw 



I. INTRODUCTION 

Theories unifying gravity with other interactions sug- 
gest the possibility of spatial and temporal variation 
of fundamental physical constants (VFC), such as the 
fine structure constant, a = e^ /fic, and the proton-to- 
electron mass ratio, /i = nip/me [l|. Search for such vari- 
ation has received considerable interest in recent years, 
and is being conducted using a wide variety of meth- 
ods [2,|3|. Some major directions of this research include 
analysis of high resolution spectroscopy of quasar absorp- 
tion systems [J-lfil , frequency comparison of atomic clocks 
over extended periods of time |7H9|, and nuclear meth- 
ods, including study of nucleosynthesis, alpha and beta 
decay, and Oklo natural reactor |10l-[la|. 

Precision molecular spectroscopy is a new and promis- 
ing direction of search for variation of fundamental con- 
stants. Molecular spectra are sensitive to both /z and a, 
and by measuring close lying levels great enhancement 
of relative variation may be observed [J, [la, [l3| • In par- 
ticular, diatomic molecules that have a near cancellation 
between hyperfine structure and rotational intervals or 
between fine structure and vibrational intervals are of 
interest in the context of such an enhancement. A num- 
ber of such molecules have been proposed, e.g. Cs2 p^ . 
CaH, MgH, CaH+_p,[23l, ClJ, IrC, HfF+, SiBr, LaS, 
LuO, and others |2]| . 

In this paper, we conduct a detailed study of one of 
the molecular candidates suggested by Flambaum and 
Kozlov [2l|, namely silicon monobromide. To this end, 
it is useful to start by briefly recapitulating some of the 
main concepts put forth by these authors. We consider a 
diatomic molecule with an electronic ground state com- 
posed of a fine structure multiplet. Taking the vibra- 
tional energy spacing of the multiplet as uie in the har- 
monic approximation, and the fine structure (spin-orbit) 
energy spacing between two multiplet states as w/, the 



energy associated with a transition between the multiplet 
states reads 

where v represents the change in the number of vibra- 
tional quanta for the transition. 

The fine structure and vibrational energies have differ- 
ent sensitivities to variations in a and fi. In particular, 
ijf is sensitive to variations in the fine structure con- 
stant, scaling as ojf ~ a^, while being almost insensitive 
to variations in fi. On the other hand, Wg is sensitive to 
variations in the proton-to-electron mass ratio, scaling as 
We ~ /i~^", while being insensitive to variations in a. It 
follows that id is sensitive to variations in both a and ^, 
with a corresponding variation for fractional variations 
in a and ^ given by 
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For a number of molecules there exist transitions hav- 
ing a near cancellation between fine structure and vi- 
brational energies, i.e., w/ « vuJe- In such cases, the 
fractional variation of a; may then be written 

where K = uif/co is an enhancement factor. Large values 
of K are suggestive of favorable cases for experimentally 
detecting a signal from variations in a or fi. As discussed 
in Rcf. [2H, however, it is also necessary to consider the 
size of the absolute shift Suj and compare this to exper- 
imental limitations on measuring w itself; one such no- 
table limitation is the natural linewidth and intensity of 
the transition. 

The diatomic molecule SiBr has a ^11^ electronic 
ground state with fine structure and vibrational spacing 



similar to about 1 cm"^ (w/ w We ~ 420 cm~^). This 
is comparable to the rotational constant B^, and thus 
Lo may be reduced further by a suitable choice of rota- 
tional levels. In this paper we examine the rovibronic 
spectrum of SiBr by employing the most accurate exper- 
imental spectroscopic data for SiBr available in the liter- 
ature, namely that of Bosser et al. [22|. Furthermore we 



perform ab initio molecular calculations with the purpose 
of i) determining the precise dependence of the spectrum 
on a, and ii) obtaining values for the natural linewidths 
of the pertinent levels. As in Ref. [2l|, we still conclude 
that dedicated measurements are required to determine 
precise values of transition frequencies and find the best 
transitions for the search of VFC; the aim of this work 
to entice experimental progress in this direction. 

At the risk of being overly prudent, we discuss our 
convention used throughout concerning units, applicable 
to the above expressions as well. We choose to work 
with atomic units {h — e — nie — 1), and thus an ex- 
pression such as 6X indicates a variation in X when ex- 
pressed in atomic units (this is not a trivial remark: for 
instance, when expressed in atomic units the speed of 
light c = l/a certainly varies with a variation of a; how- 
ever by definition the speed of light does not change if 
expressed in SI units). Throughout this paper we will 
find it useful to express energy values in the spectroscop- 
ically familiar units of cm~^; one should interpret this 
merely as a conversion from the atomic unit of energy, 
1 a.u. = 2.19474625 x 10^ cm"^ In the end we wiU only 
be concerned with variations of dimensionless quantities 
[l|, such as the ratio of two frequencies, and for these 
expressions ambiguity surrounding units is non-existent. 



II. ROVIBRONIC ENERGY LEVELS IN 
HUND'S CASE a DIATOMICS 

We consider an electronic multiplet of a diatomic 
molecule which is categorically described by Hund's case 
a |2j]. In Hund's case a, the electronic orbital angu- 
lar momentum L is strongly coupled to the internuclear 
axis (chosen to be the z-axis in a molecule- fixed frame), 
which is to say that A, the eigenvalue of L^, remains 
a good quantum number. Furthermore, the spin angu- 
lar momentum S is strongly coupled to the internuclear 
axis by way of the spin-orbit interaction, and thus S, the 
eigenvalue of Sz, also remains a good quantum number. 

Initially we neglect the spin-orbit interaction, in which 
case we may write the vibronic energies T^v of a given 
electronic multiplet in terms of conventional spectro- 
scopic constants. 
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where v is the vibrational quantum number and terms 
beyond second order in {v + i) are omitted. The con- 
stant Tg is the energy relative to the ground state multi- 
plet; as we will only be concerned with the ground state 



multiplet, we may set Tg = 0. Constants uj^ and UeXe 
represent the harmonic vibrational energy and the first 
correction due to anharmocity, respectively. 

Next we consider the effective spin-orbit interaction. 
As L is strongly coupled to the internuclear (z) axis, the 
spin-orbit interaction takes the simple form |24l | 



-f^so — ^v-^z^z 



(2) 



The spin-orbit factor A.^, here depends on the vibrational 
state and to the first order in {v + i) may be written as 



A,. 



OLAe 



1 
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(3) 



Finally, we consider the energy associated with rota- 
tion, taking the effective rotational Hamiltonian for the 
Hund's intermediate case a — 6 as in Ref. [241, 



-ffrot ~ ByN , 



(4) 



where N = J — S, and J being the total angular mo- 
mentum excluding nuclear spin. We now introduce the 
operators J = JxHJy, and similar for S^, where x and 
y correspond to the molecule-fixed axes perpendicular to 
the internuclear axis. With these operators we expand 
N^ as 



N^ 



J^ + S^-23-S 



J' 



S^ - 2Sl - 2LzSz 



[j+s- + j-s^ 



(5) 



where we have used Jz — Lz — Sz ~ with the physical 
reasoning that the molecule rotates about an axis per- 
pendicular to the internuclear axis. In the expressions 
to follow, we neglect the small w-dependence of By and 
use By = Be- For the Hund's case a limit, where S 
is assumed to be a good quantum number, the term in 
parenthesis in Eq. ([5]) may be dropped as it involves the 
raising and lowering operators S . 

We now consider the energy levels specific to a H- 
doublct, such as the ground electronic state of SiBr. The 
appropriate basis for Hund's case a is |JAf, AS), where 
AI is the eigenvalue oi Jz, Z being a space- fixed axis. 
In this basis, the doubly-degenerate ^Ili/2 state is repre- 
sented by \JM, ±1, ^pi) and the doubly-degenerate ^113/2 
state represented by \JM, ±1, ±i). The angular momen- 
tum quantum number J is necessarily a half-integer, with 
limitations J > 1/2 and J > 3/2 for the ^lli/2 and ^113/2 
states, respectively. In terms of the spectroscopic con- 
stants introduced above, the energy levels are then given 

by 



Eyj = ±^{Ae - 2Be) + (uJe T ^aAej U' + ^ 



UJeXe [V+-] + Be [ J 



1 



(6) 



where the top (bottom) sign corresponds to the ^113/2 
(^Hi/2) levels. As discussed in Ref. [2J], aAe may be 



regarded as the difference in the harmonic vibrational 
energies of the doublet levels when considered indepen- 
dently, i.e., aAe — oje — uje ; this interpretation is 
consistent with Eq. ^. 

It will be useful to separate the energy given by Eq. (|6]) 
into J-independent and J-dependent parts. 



E; 



vj 



G,„ + Fj, 



where Fj = Be [J + -^j ■ We will refer to these as the 
vibronic and rotational energies, respectively. Note that 
with this nomenclature the energy contribution T^e of 
Eq. dni) is associated with the vibronic energy, despite 
arising from the rotational Hamiltonian, Eq. (|4]). We 
further note that with this choice of separation, only the 
expression for the vibronic part depends on the particular 
doublet state (^113/2 or ^ni/2)- The total energies, Eyj 
will be referred to as the rovibronic energies. 

It is briefly noted here that certain additional terms, 
such as those associated with the spin-rotation inter- 
action, lambda-doubling, and the hyperfine interaction, 
have intentionally been neglected in this section. These 
contributions are generally small, though for large J some 
of these terms may have a sizable effect. 



III. ROVIBRONIC ENERGY LEVELS OF 
SILICON MONOBROMIDE 

The ground electronic multiplet, X ^Hr, of SiBr falls 
into the category of Hund's case a. Accurate spectro- 
scopic constants for this doublet have been experimen- 
tally determined by Bosser et al. |22l | and are presented 
in Table HI Figure [TJa) illustrates the potential energy 
curves for the ^Tli/2 and ■^113/2 states near minima, based 
on the data for isotopic species ^^Si^^Br; also displayed 
are the lowest few vibronic energy levels. Due to simi- 
larity in the magnitude of the A^ and We constants, the 



.(1/2) 



.(3/2) 



Gv level is quasi-degenerate with the G^_]^ level for 
V ~ 1,2, . . . . Figure [IJb) provides a magnification of the 
energy separation between these quasi-degenerate levels. 

Before proceeding we briefly discuss the accuracy of 
the experimental data in Table HI Explicit uncertainties 
are n ot p rovided for the constants, but indications from 
Rcf. [221 are that the data are likely to be accurate to 
~0.1 cm""'^. For the illustrative purposes of this section 
we will treat the data as exact; for deviations on the order 
of ^ 0.1 cm~^ the important qualitative features of the 
spectrum remain, and only minor modifications would be 
necessary. 

As mentioned in the Introduction, we are interested in 
transitions having large enhancement factors, namely the 
transitions between the quasi-degenerate vibronic lev- 
els. We define the small energy difference between quasi- 
degenerate vibronic levels as 



AG, = G, 



(3/2) 



(-(1/2) 



TABLE I. Spectroscopic constants for the X ^Ilr ground dou- 
blet of SiBr. Theoretical values for ^^Si^^Br are calculated 
using the relativistic Fock space coupled cluster approach, de- 
scribed in Section V. Experimental data for isotope '^^Si^^Br 
are from Ref. [22l, whereas data for isotope '^^Si^^Br are in- 
ferred using the appropriate dependence on reduced mass per 
spectroscopic constant along with the ratio of reduced masses 
of the two isotopic species, 1.0065044. All values in the table 
are in cm~^. 





^«Si'^«Br 


^«Si''«Br 


^«Si«^Br 


Constant 


theor. 


expt. 


expt. 


Ae 


419.54 


422.61 


422.61 


UJe 


424.35 


424.14 


422.77 


OJeXe 


1.32 


1.41 


1.40 


ClAa 


2.26 


1.97 


1.96 


B, 


0.1634 


0.1671 


0.1660 




CiJe — 2_Be -|- V {2u)eXe — OLAe) 



(7) 



Nuclear Separation 



FIG. 1. (Color online) Potential energy curves and the lowest 
few vibronic energy levels for the X ^lir ground doublet of 
28gj79gj. j-g^-j rpj^g lower blue and upper red curves represent 
the potential energy curves for the '^^\/2 and ^113/2 states, 
respectively, with the solid blue and dashed red horizontal 
lines illustrating the corresponding vibronic energy levels, (b) 
The energy separation between the quasi-degenerate vibronic 
levels is magnified by a factor of 100. All energy differences 
are in cm^^. 



(It is noted that AG, as defined here is not related to the 
usual spectroscopic AG1/2, AG3/2, ■ • ■ ) For the isotope 
28gj79gj. ^Yi\s reduces to 

AGi, = (-l.86-h0.85i;) cm^^ 

An interesting property of the ^^Si^^Br vibronic spectrum 
is that AG, is negative for v = 1,2 and positive for 
1) > 3. This is plainly seen in Figure [IJb), where for the 



(5/2) 

quasi-degenerate levels described by u = 1,2 the Gl,_i 

energy (dashed red line) is below the Gv energy (solid 
blue line), whereas the order is inverted for v > 3. This 
inversion arises due to the anharmocity of the potentials, 

UleXe- 

We now turn our attention to rotational energies. With 
our choice of separation for "vibronic" and "rotational" 
contributions to the total energy, the rotational energies 
are given by the same expression for both doublet states. 



p(l/2) 



^(3/2) 



I.e., Py^' = Py^' = BeiJ + ^y. We win concern 



ourselves only with single-photon transitions, from which 
the angular momentum restriction A J = 0,±1 follows. 
For A J = transitions, there is no change in rotational 
energy, and the corresponding measured transition lines 
for all J may be blended, limiting the accuracy. For this 
reason, we focus on transitions with A J = ±1. We define 
the difference in rotational energy encompassing both of 
these cases as 



AF^ 



F 



(3/2) 
J±l 



F 



(1/2) 



2Be (J H 
-2B,J 



1) 



(8) 



where AFj'" and Ai^J are restricted by J > 1/2 and 
J > 5/2, respectively. We note that AF^ is necessarily 
positive, whereas AFJ" is necessarily negative. 

The experimentally observable quantity is the energy 
difference between two rovibronic levels; we define the 
energy difference between pertinent rovibronic levels as 



^Etj 



E, 



(3/2) 



-1,J±1 



E 



(1/2) 
vJ 



AG, 



AFf. 



(9) 



To continue with our strategy of finding the largest en- 
hancement factors, we look for specific transitions in 
which 



line) is above the E\^j ' level (solid blue line), then the 
sign of ^E^j is positive; if the order of the levels is op- 
posite, the sign of AF^j is negative. The particular sign 
of AFy J is important from the viewpoint of variations of 
AFyj with respect to variations of a and /i, as discussed 
in the following section. 



a) 



b) 



■ J=9/2 



ni/2, v=i 




^n3/2,V = 


t 

1.01 

1 , 



/=9/2 

J=7/2 

-f -^=5/2 

tQ-18 J=3/2 
-^ J=5/2 

y= 1/2 

J=3/2 



FIG. 2. (Color online) Rovibronic levels of ^*Si^^Br associated 
with the ^Tlif2, v — 1 (solid blue lines) and ^113/2, v = 
(dashed red lines) quasi-degenerate vibronic levels. Displayed 
are (a) the two vibronic levels and (b) all rovibronic levels for 
up to J = 9/2. All energy differences are in cm~^. 



AF± = AG, 



AF± 



0. 



As an example, we take the v = 1 vibronic energy differ- 
ence of ^SSi^^Br, AGi = -1.01 cm^i. As AGi is nega- 
tive, we require AFj!' in Eq. (jlOp and may subsequently 
solve for J, 



Ji 



AGi 
2F, 



-1 = 2.02, 



which indicates that two appropriate choices for J are 
J = 3/2 and J = 5/2, with corresponding values of AF^j 



(10) IV. VARIATIONS OF THE ROVIBRONIC 

TRANSITION FREQUENCIES WITH RESPECT 
TO VARIATIONS OF a AND ai 

In this section we consider variation of the energy dif- 
ference ^E^j with respect to variations of a and //. The 
constants Ae and Ue are orders of magnitude larger than 
the other spectroscopic constants used to describe ^E^j; 
furthermore A,, is only sensitive to variations in a while 
We is only sensitive to variations in /i. Consequently, 



AF+3/2 = -0.18 cm" \ 
AF+5/2 = +0.16 cm" \ 



(11) 



Figure [5] displays the rovibronic spectrum arising from 



(1/2 



the G\"^> (^ni/2, 



1) and G, 



(3/2) 



('n3/2 



0) 



quasi-degenerate vibronic levels for ^Si Br. The en- 
ergy differences corresponding to AF+g ,2 ^"^^ ^^1^5/2 
are also displayed. One can resolve the absolute values 
appearing in Figure [5] with the signed values appearing 
in Eq. ([TT|) by noting that if the -E,„_j j_j_-^ (red dashed 



to a first approximation we can estimate the variation 
5 (AFy j) by variations of A^ and 



LOp 



5{^Efj)'^5A,-5uje. 

The spin-orbit constant A,, enrbodies the major rela- 
tivistic correction to the energy spectrum of the doublet 
and to the lowest order scales as a^ . Thus, if we assume 
higher-order relativistic corrections to be negligible, we 
may write 

Me=2Ae — . 



The harmonic vibrational energy We is insensitive 
to relativistic corrections, though it is proportional to 

— 1/2 

■^red 1 where M^cd is the reduced nuclear mass, and as 
such is sensitive to /x. The proton and neutron masses, 
as well as nuclear binding energies, are all proportional 
to the quantum chromo dynamics scale Aqcd (see, e.g., 
Rcfs. [2^ [23 )• It follows that the nuclear masses and, 
further, the reduced nuclear mass are also proportional 
to Aqcd- We conclude that JM,.od/^jfrcd = Smp/rup = 
Sfi/ fi, where the last equality holds for atomic units. The 
constant uje then varies with /i as 



Suje 



1 Sfj. 

2 // 



Combining the above equations yields 

r / . + \ , 5a 1 Su 



a 

„ , / Sa I Su 

= 2AA — + 7 — 

a 4 /i 



(12) 



where in the last expression we used the fact that Ae = 
uje- Evidently the transitions AiJ^j are sensitive to 
variations in the combined constant a^'^ (the term in 
parentheses is equivalent to the fractional variation 5p/ p, 
where p = ap}''^). 

As discussed in the Introduction, the variation 
5 (Ai?,y j) is dependent on our choice of unit system, 
namely atomic units. To remove dependence on the unit 
system we consider variations of dimcnsionless quanti- 
ties, such as the ratio of two transition energies. In the 
present set-up, we consider two transition energies within 
the same doublet rovibronic spectrum. Equating oji and 
UJ2 to two separate transition energies AE^j, the varia- 
tion of the dimensionless ratio ujxjuj^ is given by 



b(i>J\lijJ2) 5ijJ\ SlU2 

{UJ1/UJ2) Wl ^2 



1 

UJl 



1 



0J2 



2A„ 



5a 
a 



I5p 
47 



The sensitivity to variations in ap^''^ is maximized by 
selecting transitions which have small values of ui and 
UJ2 and which additionally differ in sign. 

Alternatively, one may measure the ratios (|a;i| -t- 
|w2|)/wr and (jwil — \u}2\) /i^r, where lOj. is some refer- 
ence energy. This is similar to what has been done with 
atomic Dysprosium (in which case wi and L02 are defined 
by different isotopic species) [23, ^M, ■ Noting that in the 
present case wi -I- W2 is sensitive to variations in ap^''^ 
whereas ui — 0^2 is not, it follows that whether or not 
|wi| ± \ijj2\ is sensitive to variations in ap^'^ depends on 
the relative signs of wi and W2- In particular, given op- 
posite signs of wi and uj2 the difference jwij — \u)2\ is sen- 
sitive to variations in ap^''^] an immediate benefit of this 
choice is that systematic effects (e.g., constant frequency 
shifts) can presumably be controlled or eliminated to a 
large extent by taking the difference. Dependence of the 
reference energy on the variation of the fundamental con- 
stants may be neglected if there is no relative enhance- 
ment (cancellation of different contributions) there. This 



is the case for the Cs hyperfinc standard and any other 
hyperfine transition (calculated in Ref. [23|) and practi- 
cally for any other transition in SiBr. 

The preceding expressions of this section are approx- 
imate relations. In the remainder of this section, we 
present more precise formulae. In particular, we include 
contributions from variations in the additional spectro- 
scopic constants WgCCe, aAe, and B^ and, further, we use 
ab initio calculations to determine the precise depen- 
dence of the variations of the constants with variations of 
a. En route to calculating the variations of the constants 
with respect to a, we obtain values for the constants 
themselves, which are presented in Table U alongside the 
experimental data. We note an impressive agreement 
between our computed constants and the experimental 
constants. In particular, the primary constants A,, and 
We agree to better than 1%. and we feel that this is in- 
dicative of the accuracy of our computed variations of 
the constants with respect to a presented below. A brief 
discussion of the computational method is provided in 
the following section. 

From our calculations we obtain the following relations 
for variations of the constants with respect to variations 
of a only 

Me 
^e 



5aAe 

OlAe 

5uJe 



5a 
2.019—, 
a 

5a 



1.927 



-4.5 X 10" 



,5a 



(13) 



and variations in constants ujeXe and Be with respect to 
variations in a give negligible contribution. 

For variations with respect to /i, we make use of ana- 
lytical formulae, using the appropriate dependence of the 
constants on the reduced mass. In particular, we have the 
following relations for variations of the constants with re- 
spect to variations of p only 

5uie 1 5p 

CJe 2 // ' 

5aAe 1 5p 

UAe 2 p ' 

5 [uJeXe] _ _5p 
UJeXe ^J■ ' 

5B, ^ _5_p^ ^^^^ 

tie M 

and no variation in Ae with respect to variations in p. 

We consider as an example the two transition energies 
LLJi = AE'^^,^ and W2 = AE^^,^ (i-e., the transitions 

identified in Fig. EJJb)). Using Eqs. (|13|14|) . we find the 
variations in these transition energies to be 



5uji = 5uj2 = (851 cm ^^ 



5oL „ ^,„5p 

h 0.247— 

a p 



The deviation of this expression from the less sophis- 
ticated expression, Eq. (H^), is small and on the order 



of the accuracy of the computations. Note, however, 
that for transitions associated with a higher vibrational 
quantum number v the deviations from Eq. (fT2|) become 
more pronounced. For example, suppose the transition 
w = Ai?^ 55 5 is found to be a convenient transition to 
probe (with Eqs. (|7|8|9|) . u = -0.01 cm"^); for this tran- 
sition, we obtain 



5uj= (764cm-M ( — 
a 



0.220— 



V. COMPUTATIONAL DETAILS OF X ^H^ 
ROVIBRONIC SPECTRUM 

Fine structure splitting is an inherently relativistic ef- 
fect; hence, all the energy calculations employed the rela- 
tivistic four-component molecular Dirac-Coulomb Hamil- 
tonian. 



i/Dc = E/^D« + E;^ 



j<j 



where 



hY){i) = ca^ ■ Pi + Ac^ + Kuc(J)- 



Here, /id is the one-electron Dirac Hamiltonian, with Vnuc 
the nuclear attraction operator for the two nuclei consid- 
ered, and takes into account the finite nucleus effect, a 
and /3 the four-dimensional Dirac matrices, and the term 
Si<7 l/''u represents the repulsive Coulomb interaction 
between electrons. 

As SiBr is an open shell system with a single valence 
electron outside the closed shell, we employed Fock space 
coupled cluster (FSCC) method with sectors (0,0) and 
(0,1) to account for electron correlation, where the closed- 
shell cation served as reference, and an electron was 
added in the (0,1) sector, with the model space composed 
of both ^ni/2 and ^113/2 molecular orbitals to obtain the 
potential energy curves for the two states of interest. 

An uncontracted aug-cc-pVTZ basis set was used for 
both atoms [30, [Sll; 37 electrons were correlated and 
virtual orbitals with energies above 35 a.u. were omit- 
ted. All the energy calculations were performed using 
the DIRAC program package [321 : ^^id the spectroscopic 
constants were obtained from the potential energy curves 
by solving the rotational-vibrational Schrodinger equa- 
tion numerically using the program VIBROT [Sj]. The 
results are presented and compared with the experimen- 
tal values in Table I. 

In order to assess the dependence of the values of in- 
terest, j4e. We, a^e: WeXe and B^,, on the fine structure 
constant, the calculations were carried out for different 
values of x = (a/ao)^ — 1, and the derivatives were ob- 
tained using numerical differentiation. 



VI. COMPUTED LINEWIDTHS OF 
ROVIBRONIC LEVELS 

The natural linewidths of the rovibronic levels place a 
fundamental limit on the accuracy to which a given tran- 
sition frequency can be measured, and, therefore, also on 
the accuracy to which any variation in a^^'^ may be mea- 
sured. In Ref. [2l|, rough estimates of the linewidths of 
the rovibronic levels were obtained for Il-doublets; here 
we provide computed values for the SiBr molecule. 

To obtain the linewidths we require dipole matrix ele- 
ments. We begin with the wave function ^^i,(g, i?) rep- 
resenting the Born-Oppenheimer molecular solutions for 
a non-rotating molecule (see, e.g., Ref. [23j). 

where R is the nuclear separation and q encapsulates 
all electronic space and spin coordinates. Here 'ip^{q,R) 
represents the electronic eigensolutions for a given nu- 
clear separation R, and 4>-yy (R) represents the subsequent 
eigensolutions for the nuclear vibrational motion. The 
solutions are assumed to be orthogonal and normalized 
over the appropriate space, i.e., 

il!*,{q,R)il;^{q,R)dq^ Sy^, 

<l)*^^,{R)(j)^y{R)dR = 5^1 y. 

We may write the dipole operator in terms of electronic 
and nuclear contributions, 

D(5,i?)=De(g)+D„(i?). 

Specifically, the electronic contribution De(g) is given by 

De(9)--E'^" 
i 

where r^ is the position vector of the i-th electron (r^ e q) 
and the summation runs over all electrons. The nuclear 
contribution D„(i?) is 

where Mi and Zi are the mass and atomic numbers of 
the nuclei i = A^B and we have assumed the coordinate 
origin to be at the center of mass with z-axis aligned with 
the internuclear axis. A dipole matrix element between 
the Born-Oppenheimer wave functions reads 

(7V|D|7«) = JJrY{q,RWy.'{R) [De(9) + D„(i?)] 
X i/jry [q, R)4>^v (R)dqdR 

c^;,,.[R){i\D\-i)cj,^y{R)dR, (15) 



where the i?-dependent matrix clement (7'|D|7) repre- 
sents the dipole matrix element for a given "clamped" 
nuclear separation R and is given by 



{R)+ / ry{q,R)^eiq)Mq,R)dq- 



(7'|D|7) = Sy^Br, 



We are interested in dipole matrix elements between 
the vibronic states of a (Hund's case a) electronic mul- 
tiplet. The relevant vibrational wave functions (f>-yy{R) 
are largely independent of the particular multiplet level 
(i.e., (j)jviR) = (f'viR))- Furthermore, the wave func- 
tions (jiviR) are significant in magnitude only within a 
small region about the equilibrium nuclear separation 
R = Re (note that the Born-Oppenheimer approxima- 
tion is rooted in the assumption that (j)v{R) varies more 
rapidly with R than 4'-({q, R))- As such, we may expand 
(7'|D|7) in the integrand of Eq. (jTSl) about R^; explicitly 
to the first order in {R — R^) this is 



(7'|D|7) - (7'|D|7) 



d(7'|D|7) 



dR 



{R — Re 



Subsequent evaluation of the integral gives 

, rf(7'|D|7) 



(7V|D|7z;)=<5„,.(7'|D|7) 



R, dR 

xRe\ —^ {^/v6v' + l,v + y^Sv'-l,v) , 
y UJe ^ ' 

(16) 

where we have assumed the vibrational wave functions 
(j}viR) to be harmonic oscillator eigenfunctions. 

The contribution to the natural linewidth for a given 
decay channel jv — >■ j'v' is given by 



r(7w 



7V)^^%^|(7V|D|7^;)f 



(17) 



where w^ 



is the energy difference between the ini- 



tial and final state, and we have summed over final rota- 
tional states. We neglect decay channels within a given 
vibronic level and between quasi-degenerate vibronic lev- 
els for which the energy difference is small. 

We begin by considering vibrational decay. Using 
the MOLPRO computational package [SJ], we calcu- 
lated the diagonal matrix elements (^ni/2|Dpni/2) and 
(^113/2 |Dpn3/2) at multiple nuclear separation distances 
R within the vicinity of Re- With numerical differentia- 
tion and taking the experimental ratio y^Be/iOe = 0.020, 
we obtain the results 



('ni/2,«- 
('n3/2,« 



l|D|2n3/2,z;) 



0.12Vwa.u. (18) 



The decay channel ^113/2, w — > '^^1/2^'" is forbidden in 
the non-relativistic limit, though it is opened up by spin- 
orbit mixing of the '^lli/2 state with the excited electronic 



^Ei/2 state. Due to its purely rclativistic origin, the cor- 
responding dipole matrix element proves very difficult to 
obtain by computational methods to any degree of relia- 
bility. To circumvent the need for a direct computational 
value, we write the dipole matrix element as in Ref. [2l|, 



^n^ 



3/2 



Dpni/2)|=C| 



'n3/2|Dp 



Sl/2 



where ^ is the small parameter quantifying the spin-orbit 
mixing. We find from computation | (^n3/2|DpI]i/2) | ~ 
0.1 a.u., indicating that |(^n3/2|Dpni/2)| is apprecia- 
bly smaller than the dipole matrix elements in Eqs. (|18p . 
Thus we conclude that neglect of the decay channel 
^113/2, w — > ^ni/2i^ is acceptable as its contribution to 
linewidths will be overshadowed by the contribution from 
the vibrational decay. 

As a specific example, we consider the natural 



linewidth of the ■^fl 
Eqs. dni, we find 



1/2, 



1 vibronic level. With 



r(2n 



l/2i 



1) = 5.3 X 10"" a.u. = 0.35 Hz. 



We note that this is an order of magnitude larger than 
the estimate given in Ref. |21j . 



VII. CONCLUSION 



Here we have extended upon the Flambaum and Ko- 
zlov's work 21| by considering properties of silicon mono- 
bromide that make it a prospective candidate for detect- 
ing variations in the fine-structure constant a and the 
proton-to-clectron mass ratio fi (in particular, variations 
in the combined constant afj,^''^). We have examined the 
rovibronic spectrum by employing the most accurate ex- 
perimental data available in the literature, namely that 
of Bosser et al. [23l . Furthermore, we present results of 
ab initio calculations for the precise dependence of the 
spectroscopic constants on variations in a. We addition- 
ally present calculated values for the natural linewidths 
of the rovibronic levels which place a fundamental limit 
on the accuracy to which variations in afi^' may be de- 
termined. 



As in Ref. [2l|, we emphasize that dedicated measure- 
ments arc necessary to find precise values for the transi- 
tion frequencies and determine the best transitions. It is 
our hope that this work entices experimental progress in 
this direction. 
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